Mon. Not. R. Astron. Soc. 000,[TV?? (0000) Printed 12 December 2012 (MN KTeX style file v2.2) 



The accretion flow to the intermittent accreting ms pulsar, HETE 
J1900.1-2455, as observed by XMM-Newton and RXTE 

"(N A. Papitto 1 *, A. D'Ai 2 , T. Di Salvo 2 , E. Egron 3 , E. Bozzo 4 L. Burderi 3 , R. Iaria 2 , 
O ; A. Riggio 3 , M. T. Menna 5 

^s^j . 1 Institut de Ciencies de VEspai (IEEC-CSIC) Campus UAB, Fac. de Ciencies, Torre C5, parell, 2a planta 08193 Barcelona, Spain 
' 2 Dipartimento di Fisica, Universitii di Palermo, via Archirafi 36, 90123 Palermo, Italy 
^ ' 3 Dipartimento di Fisica, Universita di Cagliari, SP Monserrato-Sestu, Km 0. 7, 09042 Monserrato, Italy 
JIT , 4 ISDC Science Data Center for Astrophysics of the University of Geneva, chemin dEcogia 16, 1290 Versoix, Switzerland 
5 1NAF Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monteporzio Catone, Italy 



12 December 2012 



ABSTRACT 

We present a study of the accretion flow to the intermittent accreting millisecond pulsar, 
HETE J1900. 1-2455, based on observations performed simultaneously by XMM-Newton and 
RXTE. The 0.33-50 keV energy spectrum is described by the sum of a hard Comptonized 
component originated in an optically thin r ~ 1 corona, a soft thermal kT lu ~ 0.2 keV com- 
ponent interpreted as accretion disc emission, and of disc reflection of the hard component. 
Two emission features are detected at energies of 0.98(1) and 6.58(7) keV, respectively. The 
latter is identified as Ka transition of Fe XXIII-XXV. A simultaneous detection in EPIC-pn, 
EPIC-MOS2, and RGS spectra favours an astrophysical origin also for the former, which has 
an energy compatible with Fe-La and helium-like Ne-Ka transitions. Broadness of the two 
features, cr/E ~ 0.1, suggests a common origin, resulting from reflection in an accretion disc 
with inclination of (30I3) , and extending down to R m = 251}^ gravitational radii from the 
compact object. However, the strength of the feature at lower energy measured by EPIC-pn 
cannot be entirely reconciled with the amplitude of the Fe Ka line, hampering the possibility 
of describing it in terms of a broad-band reflection model, and preventing a firm identification. 
Pulsations at the known 377.3 Hz spin frequency could not be detected, with an upper limit 
of 0.4 per cent at 3-a confidence level on the pulsed fractional amplitude. 

We interpret the value of the inner disc radius estimated from spectral modelling and the 
lack of significant detection of coherent X-ray pulsations as an indication of a disc accretion 
flow truncated by some mechanism connected to the overall evolution of the accretion disc, 
rather than by the neutron star magnetic field. This is compatible with the extremely close 
similarity of spectral and temporal properties of this source with respect to other, non pulsing 
atoll sources in the hard state. 

Key words: line: identification - line: profiles - pulsars: individual HETE J1900. 1-2455 - 
X-rays: binaries - stars: neutron - 



1 INTRODUCTION 

Discovery of millisecond coherent pulsations from a number 
of neutron stars (NS) in lo w mass X-ray binaries (LMXB, 
Wijnands & van der Klisl[l9 98) proved that accretion of disc an- 
gular momentum is an efficient mechanism to spin up a NS to 
ms spin periods. X-ray pulsations from an accreting NS are ob- 
served when the NS magnetic field is strong enough to control at 
least partly the mass flow in the NS surroundings, channelling it 
towards the magnetic poles. So far, however, coherent ms signals 
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have been detected only from a subset of LMXB; o f the ~ 60 
systems confirmed to host a NS dRitter & Kolbll2003l update RK- 
cat7.18, 2012), coherent pulsations with a period in the ms range 
were observed from only 14 sources, commonly known as a ccret- 
ing millisecond pulsars (AMSP; see IPatruno & Wattsll2012l for a 
recent review). Whether a magnetosphere is able to control mass 
flow around the NS depends on the balance between the torques 
exerted on the accretion flow by the magnetic field and by disc 
matter viscosity. Observation of ms pulsations, thermonuclear X- 
ray bursts, and quasi periodical oscillations, strongly points to- 
ward the magnetic field of most NS in LMXB to lie in the range 
10 8 — 10 10 G. Lack of a detection of pulsations from majority 
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of LMXB, with present-day observatories, is thus most easily in- 
terpreted in terms of a larger average mass accretion rate; this en- 
hances diamagnetic screening of NS magnetic field by the accretion 
flow, and burial of the field under the NS surface (Cummi ng et al.l 
1200 lh . AMSP are indeed relatively faint X-ray transients. Most of 
them show few-weeks to few-months long outbursts, separated by 
years-long intervals spent in a quiescent state. Their peak outburst 
X-ray luminosity never exceeds ~ 10 37 erg s _1 , and the long- 
term average mass accretion ra te is usually ;$ 10 -3 times the Ed- 
dington rate (Gallowav 2006). However, lack of detec ted pulsa- 
tions from very faint accreting NS (e.g., lPatrunoll2010l) indicates 
how such a picture is not entirely satisfactory; alternative scenarios 
proposed to explain the paucity of ms pulsars in LMXBs point to 
puls e decoherence induced by scattering in a hot cloud around the 
NS (Titarchuk et al. 2002), m agneto-hydrodynamical instabilities 
dKulkarni & Romanovj l20oj ; iRomanova etai] [2008), alignment 
betwe en the magnetic field and spin axes of the NS dLamb etai] 
2009), a nd a low pulse amplitude ca used by gravitational light 
bending dWood et al.lll988l ; l6zelll2009h . 

A source of crucial importance to investigate differences be- 
tween pulsating and ordinary LMXBs is HETE J1900. 1-2455. Dis- 
covered in 2005 by High Energy Transient Explorer 2 (HETE-2, 
ISuzuki et alj|2007l) . Rossi X-ray Timing Explorer (RXTE) d etected 
pulsations at a frequency of 377.3 Hz dKaaret et alj |2006). Sub- 
sequent observations showed how pulsed fractional amplitude de- 
creased from ~ 4.5 per cent to below the sensitivity level, ,$ 1 
per cent. This so urce was then the firs t AMSP to show inter- 
mi ttent pulsations dGallowav et alj|2007h . later j oined by Aql X- 
1 dCaseiraetalJl2008l) and SAX J1748.9-2021 dAltamirano etai] 
2008). The magnetic burial scenario is particularly appealing to ex- 
plain the behaviour observed from HETE J1900. 1-2455, which is 
the only AMSP to have shown an outburst lasting more than seven 
years; in fact, disappearance of pulsations after two months of ac- 
cretion at a rate of few per cent the Eddin gton rate is com patible 
with the timescale of magnetic screening (Cumming 2008). After 
their first disappearance a few months after the outburst onset, pul- 
sations from this source were sporadically observe d to re-appear at 
a ver y low fractiona l amplitude, ~ 0.5 per cent dGallowav et al.l 
2008; |Patrunoll2012l) . This suggests how the magnetic field can 
temporarily re-emerge to channel accretion flow towards NS poles, 
or that emission is always pulsed, but most of the time below de- 
tectability threshold of current observations. 

A powerful probe of the accretion flow close to a compact 
object is the reflection component sometimes detected in the X-ray 
energy spectra of these sources. Such component arises from partial 
back-scattering and reprocessing of ha rd ionising X-rays il luminat- 
ing the disc (see, e.g., the review by Fabian & Ross 2010). As the 
disc matter moves in high-velocity Keplerian orbits in the gravita- 
tional well of the compact object, the shape of reflect ed spectrum 
is mo dified by Doppler and gravitational redshifts CFab ian et al.l 
1989), and conveys crucial information about the flow geometry. 
In particular, the most prominent reflection feature is an emis- 
sion line at the energy of Ka transition of iron; such feature is 
mainly due to fluorescence of mildly ionised to neutral iron (~ 
6.4 keV), or to recombination of helium-like ions of iron (6.67- 
6.7 and 6.97 keV for FeXXV and XXVI, respectively; see, e.g. 
iKallman et alj2004h . XMM-Newton and Suzaku observations of the 
AMSP, SAX J1808.4-3658, have shown how the broad shape of 
an iron line could be used to probe the transition region between 
disc an d magnetospheric flow around a quickly rotating accreting 
pulsar JPapitto et al]|2009l;ICackett et al.ll2009l ; |Patruno et alj2009l ; 
I Wilkinson et al.l201lh . In particular, the inner disc radius of the re- 



flecting layer was estimated to lie between 12 and 26 km from an 
assumed 1.4 Mq NS, compatible with the theoretical expectations 
for a NS with a few x 10 8 G magnetic field, accreting at the rate 
indicated by the observed X-ray flux. 

Here, we present the analysis of an XMM-Newton observation 
of HETE 11900.1-2455 performed in September 2011, aimed at 
determining the properties of the accretion flow to the NS on the 
basis of the shape of the disc reflection component and of its time 
variability properties. To obtain a coverage at higher energies, we 
also take advantage of a simultaneous observation performed by 
RXTE. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 XMM-Newton 

XMM-Newton d-Tansen et alj|200lh observed HETE J1900.1-2455 
for 65.3 ks starting on 2011 September 19, at 15:31 UTC (Obs.Id. 
0671880101). Data were processed with the latest version of the 
XMM-Newton Science Analysis Software released at the moment 
of performing the anal ysis (v. 12.0.0). 

EPIC-pn camera ( Strii der et al 1l200lT) was operated in Timing 
Mode to achieve a temporal resolution of 29.52 \is. We discarded 
last 5 ks of data from the analysis, since a marked increase of single 
pattern events in the 10-12 keV band indicated contamination by 
flaring particle background. The 0.6-11 keV light curve, corrected 
with epiclccorr, is shown in Fig. Q] A burst clearly appears 
above an average count-rate of 192.5 s -1 (see the inset of Fig. [T] 
where the burst light curve is plotted at a magnified scale). When 
analysing persistent emission, we discarded 5 s before and 400 s af- 
ter the burst onset. In Timing observing mode imaging information 
is preserved in one dimension only (RAWX), while on the other 
(RAWY) data are collapsed into a single row to allow a faster read- 
out. To extract a spectrum of the source emission we first consid- 
ered 0.6-11 keV photons falling in a rectangular region covering 
RAWX = 29 - 46, discarding events flagged as bad (FLAG=0) 
and retaining only single and double pattern events (PATTERNS). 
We used the task epatplot to evaluate the observed fraction of 
single and double events with respect to expectations. A devia- 
tion up to 2% with respect to the model was found in the spec- 
trum extracted from a region covering also the brightest RAWX 
columns. We could obtain a fraction of single and double events 
compatible with the model only by excising the two brightest ones 
(RAWX=37-38); in order to avoid any slight spectral distortion of 
the continuum due to photon pile-up we considered the spectrum 
extracted from this reduced area. Spectra were re-binned with the 
task specgroup, to not over-sample the instrument energy reso- 
lution by more than a factor of three, and to have at least 25 counts 
in each energy bin. The ancillary response file (ARF) was produced 
by subtracting the ARF obtained from the excluded region to the 
ARF calculated without excluding central columns. To perform a 
temporal analysis we retained events with all patterns and falling 
also in the brighter RAWX columns of the CCD. 

Of the two EPIC-MOS cameras dTurner et alj|200ll) , MOS2 
was operated in timing mode to reduce pile-up, achieving a time 
resolution of 1 .75 ms, while MOS 1 could not observe in such mode 
and was kept switched off to allocate more telemetry bandwidth 
to the other CCDs. MOS2 observed an average count rate of 55.8 
s ; we considered single pattern (PATTERN=0) photons falling 
in RAWX = 276 - 303 and 306 - 335 to extract a spectrum, excising 
the central CCD columns to limit pile-up, and removing bad events 
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Figure 1. 0.6-1 1 keV light curve observed EPIC-pn. Inset shows the burst 
light curve with the time axis shown in seconds since the burst onset. 



(FLAG=0). Background was extracted from a box-shaped, 3000 x 
2000 pixel wide region free of sources, located in one of the outer 
CCD operated in imaging mode. We used similar re-binning param- 
eters than those applied to produce EPIC-pn spectrum. We could 
not consider MOS2 data to perform a temporal analysis aimed at 
detecting a coherent signal, as the Nyquist frequency associated to 
its time resolution, 285.7 Hz, is lower than the spin frequency of 

the sorce (377.3 Hz). 

Reflection grating spectrometer (den Herder et al .1 .2001) oper- 
ated in standard spectroscopy mode. We considered only first-order 
spectra, re-binned in order to have at least 25 counts per channel. 
The average persistent count rate was 3.9 and 4.8 c s _1 for RGS1 
and RGS2, respectively. 



2.2 Rossi X-ray Timing Explorer 

Simultaneously to XMM-Newton observation, HETE J1900. 1-2455 
was also observed for 3 ks by Rossi X-ray Timing Explorer (RXTE), 
starting on 20 September at 08:50 (ObsId.:96030-01-34-00). We 
extracted a spectru m from Proportional Counter Array data (PC A 
Jaho da et al ] |2006h . considering data taken by the top layer of 
PCU2 only, and producing a response matrix with the latest cali- 
bration issued (Shaposhnikov et al. 20093)- Following their recom- 
mendations we added a systematic error of 0.5% to each spectral 
bin. A spectrum from data taken by High Energy X-ray Timing 
Experiment (HEXTE lRothschild et all 1998) was extracted consid- 
ering cluster A, as it was the only in a on-source position during the 
observation. Data observed by cluster B, which was instead perma- 
nently in a off-source position, were used to estimate background. 



3 X-RAY SPECTRUM OF PERSISTENT EMISSION 
3.1 The EPIC-pn 0.6-11 keV spectrum 

The X-ray spectrum of AMSP is dominated by Comptonization 
of soft (kT 30 ft ~ 0.5-1 keV) photons in a hot (fcT f > 20 - 
40 keV), optically thin (r ~ 1) medium JPoutanenl 120061) . 



When observations at energies below ~ 2 keV are avail- 
able, one or two thermal components are detected and com- 
monly interpreted as originating f r om th e accretion disc and/or 
the NS surface jGierlinski et al.| 2002 ; Gierlinski & Poutanenl 
120051 ; IPapitto et alj|2009l ; IPatruno et al.ll2009l : IPapitto et alj|20im 
A broad iron Kg emission feature i s often detected as we ll 
( IPapitto et alj|2009l ; ICackett et al.ll2009l |201(]| : IPapitto et allboiol) . 
and to gether with a Compton hump peaking between 20 and 40 
keV ( Ibragimov & Poutane rJ2009l ; |Papitto et al.l20ld ; lKaiava et all 
1 20 111 ; llbragimov et al. 201 ll) is generally interpreted as a finger- 



print of disc reflection of the hard X-ray emission generated close 
to the compact object. 

We first modelled the 0.6-11 keV continuum spectrum ob- 
served by EPIC-pn (see the top panel of Fig. |2) using a thermal 
Comptonized component and a black-body, and describing pho- 
toelectric absorption using an improved version of the Tiibingen- 
Boulder model (TbNew in XSpec terminology, Wilms, Juett, 
Schulz, Nowak, 2011, in preparatiorjf]). Abundances and photo- 
electric cros s -sectio ns were set following Wil ms et al. I (Hi) and 
IVerner et all £996), respe ctively. We modelled Comptonization 
with the nthcomp model (Zdziarski et al. 1996; Zvcki et al. 1999, 
see blue line in panel (a) of Fig. O, which describes the up-scattered 
spectrum in terms of a power-law of index r ~ 2, between a 
low-energy rollover at the seed black-body photon temperature, 
fcT so ft ~ 0.45 keV, and a high-energy cut-off at the electron tem- 
perature of the Comptonizing medium; since this is left uncon- 
strained even by data obtained by RXTE at higher energies (see 
Sec. [33}, we held it fixed to kT e = 50 keV. The thermal soft 
component has a temperature of ~ 0.2 keV and a normalisation 
corresponding to a radius of £ 30 ri s km, where ris is the dis- 
tanc e to the source in uni ts of 5 kpc dKawai & Suzukill2005l ; see 
also iGallowav et al.l 120081 who derived a compatible estimate of 
4.7 ± 0.6 kpc). Since this size exceeds the radius of a NS, and con- 
sidering how the temperature is typical of accretion disks around 
this kind of sources (see, e.g., IPapitto et al.l l2009h . we replaced 
the black-body component with a multicolour accretion disc model 
(diskbb in XSpec; see red line in panel (a) of Fig. |2j- Evident 
residuals at ~ 37 and ~ 5 a above the continuum model appear 
around Ex = 0.98(1) and E 2 = 6.58(7) keV. We modelled them 
with Gaussian emission lines with a width of <7i = 0.12(2) and 
a 2 = 0.7(1) keV, respectively (see Fig. [3] showing the residuals 
obtained when these two lines are removed from the model). The 
description of the spectrum is also improved by an absorption line 
centred at E3 = 1.79(2) keV, probably due to a mis-calibration of 
the instrumental Si edge (see below). The chi-squared we obtained 
after the addition of these features is 1.40 over 169 degrees of free- 
dom; residuals are shown in panel (b) of Fig. [2] and the values of 
the parameters obtained are listed in the column labelled as Model 
A in Table[D 

Given the presence of small residuals at low ^ 2 keV ener- 
gies, and of a weak absorption feature at an energy close to the 
instrumental Si edge (1.84 keV), we also modelled the spectrum 
obtained after correcting the data with the task epf ast. This tool 
is designed to reduce the possible effect of charge transfer inef- 
ficiency (CTI), and/or X-ray loading (XRL), on the spectra ob- 
served by the EPIC-pn in fast readout modes. It was calibrated 
to reduce residuals with respect to phenomenological models in 
the 1.5-3 keV, where the largest deviations are usually found (see 



1 http://heasarc.gsfc.nasa. gov/ docs/xte/pca/doc 
|/rmf /pcarmf-11 . 7/| 



2 |http : / /pulsar . sternwarte ■ uni-erlangen . de/wilms/| 
research/tbabs/ 
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Figure 2. 0.6-1 1 keV spectrum observed by EPIC-pn during the persistent 
emission. The Comptonized component (blue), the accretion disc multi- 
temperature black-body (red) and the Gaussian features at 0.98 and 6.58 
keV (magenta and light blue, respectively), used to model the spectrum ac- 
cording to Model A (see Table [T} are also plotted (panel a). Residuals of 
EPIC-pn data with respect to the best fitting Model A and B (see text and 
caption of Table[TJ are plotted in panel (b) and (c), respectively. Residuals 
of the spectrum simultaneously observed by EPIC-MOS2 with respect to 
Model A (to which a couple of Gaussian absorption feature at ~ 2 keV 
were added) are plotted in panel (d). 



the XMM-Newton calibration technical note, XMM-SOC-CAL- 
TN-0083, available at http : / / xmm . esac . esa . int /] ). How- 
ever, the corrected spectrum resulted slightly more noisy than the 
unmodified one, and a further emission feature centred at 2.23(1) 
keV, close to the instrumental Au edge, should be added to the 
model; still the reduced chi squared of the fit (1.44 over 167 d.o.f.) 
was larger than the one obtained without correction, and for this 
reason we considered only the spectrum obtained without epf ast 
correction. We note how the width and the normalisation of the 
feature at higher energy is left unchanged by epf ast, while its 
centroid energy slightly increases with respect to the unmodified 
data, 6.66(7) keV. An increase of t he iron line energy in epf ast 
corrected data was already noted bv lWalton et ail (20 12), who com- 
pared XMM-Newton and BeppoSAX spectra of XTE J 1650-500, 
concluding how the unmodified EPIC-pn data possibly provided 
a better approximation of the true energy scale. 

Reflection of the hard X-ray continuum above the accretion 
disc is the most convincing intetpretation of the broad feature ap- 
pearing at 6.58(7) keV, an energy compat ible with Kg transit ion 
of iron. We used the model ref lionx teoss & Fabianll2005l) to 
describe reflection arising from an optically thick ionised slab, illu- 
minated by a r = 2 power law spectrum. We convolved the reflec- 
tion component with an exponential low-energy cutoff (expabs 
in XSpec), at an energy equal to the temperature of the soft seed 
photons of the nthcomp component, kTsoft, to take into account 
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Figure 3. Residuals of the 0.6-1 1 keV EPIC-pn spectrum with respect to 
Model A (see text and Table[TJ, when the Gaussian emission features cen- 
tred at 0.98 and 6.58 keV are removed from the model. 



the low energy rollover of the illuminating Comptonized spectrum. 
The reduced chi-squared obtained adding the reflection component 
to Model A, and removing the Gaussian feature at the energy of 
the iron Ktv transition, is 1.45 (250 over 172 d.o.f.). The ionisation 
parameter lies in the range log £ =2.7-2.8, compatible with the ob- 
serve d centroid energy of the iron line (see Table 1 of iGarcfa et al.l 
1201 ll) . The addition of a Gaussian emission line centred at ~ 1 keV 
was still requested by the data, as the reflection model could not de- 
scribe self-consistently the soft excess; removing the line from the 
model yielded in fact to a reduced chi-squared of ~ 4.5, due to 
evident residuals at those energies. At the same time, the introduc- 
tion of the reflection component made the soft thermal component 
not significantly detected anymore. However, the EPIC-pn response 
at energies ,$ 2 keV is clearly dominated by the Gaussian-shaped 
soft excess, found at a lower energy (0.78(5) keV) and with a larger 
equivalent width (153 eV) than the one obtained with model A, in 
which disc reflection was not self-consistently modelled. The non 
detection of the soft thermal component is then most probably due 
to correlation with the emission feature, whose energy is evaluated 
by this model to lie very close to the low energy end of the EPIC- 
pn energy band (0.6 keV). This is also confirmed by the detection 
of the disc thermal component in the spectrum obtained adding the 
RGS data-set, even when a reflection component is included in the 
model (see Sec. l3.3t . 

To investigate if the broadness of the emission features de- 
tected in the spectrum (a/E ~ 0.1) could be explained in terms of 
relativistic effects expected to develop in the inner parts of the ac- 
cretion disc 2 _w£_cc^rvolyed the reflection model with the rdblur 
kernel dFabian et al.ll989h . This model describes the relativistic ef- 
fects due to the motion of plasma in a Keplerian accretion disc im- 
mersed in the gravitational well of the compact object, in terms 
of the inner and the outer radius of the disc, Ri n and R ou t (in 
units of the gravitational radius, R g — GM/c 2 , where M is the 
mass of the compact object), of the index of the assumed power- 
law dependence of the disc emissivity on the distance from the 
NS, f3, and of the system inclination, i, respectively, by assuming 
the Schwarzschild metric. We held the value of the outer radius of 
the disc fixed to 1000 i? g as spectral fitting could not constrain it. 
Convolving the reflection model with such a component, the chi- 
squared of the fit decreased by A^ 2 = —20.9, for the addition of 
three free parameters. In order to check the significance of such an 
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Table 1. Best fitting parameters of the X-ray spectrum of HETE J1900.1-2455. Columns labelled as Model A (EPn), 
TbNew (diskbb+nthcomp+gaui+gau2+gau3) , and Model B (EPn), TbNew (nthcomp+expabs*rdblur*ref lionx+gau2+gau3) , re- 
spectively, refer to the 0.6-11 keV spectrum observed by the EPIC-pn. Column labelled as Model A (RGS) refers to the modelling of the RGS1/2 spectra 
in the 0.33-2.07 keV. Column labelled as Model B (XMM-XTE), TbAbs*edge (nthcomp+expabs*rdblur*ref lionx+gau2+ gau3+gau4), 
reports the parameters obtained by modelling of the spectra obtained simultaneously by RGS1/2 (0.33-2.07 keV), EPIC-pn (0.6-11 keV), PCA (3-30 keV) 
and HEXTE (20-50 keV). Quoted fluxes are unabsorbed and evaluated in the 0.5-10 keV band for Model A (EPn) and B (EPn), in the 0.5-2 keV for Model 
A (RGS), and in the 0.5-50 keV band for Model B (XMM-XTE). The inner disc radius evaluated from the normalisation of the diskbb component, 7?disc> 
is given for an inclination of 30° . Uncertainties are evaluated at a 90 per cent confidence level. Parameters held fixed in the fitting are reported without an 
uncertainty 

Model A (EPn) Model B (EPn) Model A (RGS) Model B (XMM-XTE) 



N H (10 21 cm 2 ) 
Ao 




1 

^-o.e 
1 


0.6(2) 
1 


1.02(5) 
2.11 ±0.15 


1.05(3) 
2.27 ±0.11 


fcT in (keV) 
-Rdisc (ds km) 




0.21(1) 

3i-i±|;l 




n 7 , + 0.03 
"•"-0.02 

29.5i 6 5 ;I 


0.202(7) 


r 

fc^soft 
fcT e 

N c (10" 2 cm" 2 s~ 




1.99(1) 
0.45(1) 
50 
6.7(2) 


1.979(7) 
0.458(4) 
50 
6.0(1) 


1.99 
051 +o.06 

u - J1 -0.04 

50 

6 51+ - 5 


1.975(4) 
0.446(5) 
50 
6.32(1) 


E\ (keV) 
o\ (KeV) 

/Yi (10 -3 cm" 2 s~ 
EWi (eV) 




0.98(1) 

U. 12(1 ) 
4.1±0.8 
42 


n 7 o+0.04 
U - /S -0.06 
U.24(J ) 

14+ 5 
1 -3 

153 


0.99(1) 
U.li(l) 

o o + l.O 

J - 8 -0.8 
38 


0.94(1) 

A 1 Q ( 1 ^ 

U. 13(1 ; 
4.3(5) 
40 


E 2 (keV) 
o-2 (keV) 

N 2 (1CT 3 cm~ 2 s~ 
EW 2 (eV) 




6.58(7) 
0.7(1) 
0.6(1) 
135 








Ri a (GM/c 2 ) 
Rout (GM/c 2 ) 

P 

i(°) 






23±™ 
1000 

27±3 




251- 
1000 
-3 7+ 1 - 1 

■'■ -2.5 

30.6l^ 


log(0 

R = f re fl/F nt h c 






276 +o.06 
0.10(1) 




0.086(7) 


E 3 (keV) 
a-3 (keV) 

iV 3 (10" 4 cm- 2 s~ 




1.79(2) 


-1.1(3) 


1.82(2) 


-1.2(3) 




1.80(2) 


-1.2(1) 


F (10" 10 erg cm" 2 

Xr (d-O.f.) 


s" 1 ) 


7.3(2) 
1.40(169) 


6.97(8) 
1.36(169) 


2.46(2) 
1.16(4201) 


13.59(8) 
1.25(4460) 



improvement we used the method of posterior predictive p-values 
dProtassov et al.l2002l) . We created 1000 fake spectra using the best 
fit model without disc smearing, and fitted them with a model both 
including and excluding the blurring component. The probability 
that the chi-squared decrease observed in real data after introduc- 
ing disk blurring is due to chance, is equal to the ratio between the 
number of cases in which a chi-squared decrease equal or larger 
than the one observed in real data was observed in fake spectra, 
and the number of trials. A chi-squared decrease equal or larger 
than that observed in real spectra was never obtained. We then esti- 
mated the probability that the improvement in spectrum description 
obtained by adding the disc smearing component to the model were 
due to chance as less than 10 -3 , concluding how this component 
was significant at a confidence level larger than 3-cr.The best-fitting 
parameters we obtained with this model, labelled as B, are given in 



the second column of Table[T] and the residuals with respect to the 
model are shown in panel (c) of Fig. [2] 



3.2 Soft excess, MOS2 and RGS spectra 

In the previous section, we described phenomenologically the 
strong excess appearing at ~ 1 keV in the EPIC-pn spectrum, 
by using a broad Gaussian emission line. Alternative modelling 
in terms of a bremsstrahlung continuum (bremss in XSpec) and 
an emission spectrum from hot diffuse gas (mekal), did not gave 
successful results in describing the excess, yielding a reduced chi- 
square of \r = 5.7 and 3.0, over 170 d.o.f., respectively. Several 
authors reported a similar excess at low energies in observations 
perfo r med by the EPIC-p n in fast modes (see, e.g. jHiemstra et all 
l201lUWaltonetalJl2012h . The amplitude of this excess seems to 
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be correlated with the column density of the interstellar absorption 
towards the considered source, as it is expected in the case of an in- 
sufficient redistribution calibration (Guainazzi et al. 2012, XMM- 
SOC-CAL-TN-08S- However, the column density of the inter- 
stellar absorption towards HETE J1900. 1-2455 (~ 10 21 cm" 2 ) is 
lower than that affecting sources for which no significant excess 
was found (e.g., SAX J2103.5+4545, which has a jV H ~ 6 X 10 21 
cm -2 ), and thus no calibration induced excess should be expected 
if a correlation with the column density holds. Therefore, despite 
uncertainties in the instrument calibration may be responsible for 
at least part of this excess, we also explored plausible physical in- 
terpretations. 

The centroid energy obtained with model A, 0.96(1) keV, is 
consistent with Ka emission of NelX-X (0.92-1.02 keV), and 
the Fe-La complex, which for log£ ranging between 2 and 3 is 
made by transitions with energy between 0.8 and 1 .2 keV dKallmad 
1 19951) . Such transitions are produced by disc reflection, and are al- 
ready included in the model reflionx; in particular the promi- 
nence of the Fe-La line is expected to grow when the re flecting slab 
is iron-rich jRoss & F abian 2005| iFabian et al. 1 12009). Removing 
the Gaussian feature and letting the iron abundance of the reflection 
component of Model B free to vary, we obtained a best-fit abun- 
dance of 2.00(5) with respect to Solar values, while other parame- 
ters of the reflecting medium were compatible with those listed in 
TablefflOa < 4.6, R in = 27± 3 R g , i = 29(1)°, log£ = 2.70(2), 
R = F le a/F n thc = 0.10(1); the reflection fraction R is evalu- 
ated as the ratio between the unabsorbed fluxes in the reflected and 
Comptonized component, in the 0.5-10 keV band). However, the 
reduced chi-square of the fit, 2.34 over 169 degrees of freedom, is 
much larger with respect to Model B as residuals appeared at ener- 
gies below 2 keV even after the addition of a diskbb component 
to the model. We also tried to fit the spectrum obtained discarding 
energies above 5 keV; a successful modelling is found with sim- 
ilar parameters than before, but this time with a larger reflection 
fraction, R = 0.27(2). The large flux in the ~ 1 keV feature with 
respect to the iron line is then the most probable reason for the 
non-achievement of a simultaneous modelling of the hard and soft 
part of EPIC-pn spectrum in terms of the same reflecting environ- 
ment. At the same time, we could not model successfully the whole 
EPIC-pn energy band even by using two reflection components. 

We then explored if part of the excess could be interpreted to 
an overabundance of Ne in the same reflector originating the Fe 
Ka line. Unfortunately, abundances of metals other than Fe cannot 
be adjusted in reflionx. We then added a line blurred by the 
same rdblur kernel used to convolve the reflection component, 
to the previous model. A line centred at 0.95(1) keV improved the 
modelling, but the chi-squared we obtained, 1.45 over 167 d.o.f, 
was still larger than that obtained simply modelling the feature with 
a Gaussian (1.36 over 169 d.o.f., see Model B listed in Table[T](, 
indicating how such a modelling was not completely satisfactory. 

In order to verify the origin of such an excess, we considered 
spectra observed by other instruments on-board XMM-Newton. The 
MOS2 spectrum is much more noisy with respect to EPIC-pn, prob- 
ably because of a more uncertain energy calibration when operated 
in timing mode; two additional absorption features were added in 
the Au K-edge region around 2 keV, but still the final reduced chi- 
square was far from being acceptable (2.62 over 166 d.o.f; see bot- 
tom panel of Fig.|2]for residuals with respect to Model A). Never- 
theless, we detected a Gaussian-shaped feature with parameters en- 
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Figure 4. Residuals of the 0.33-2.07 keV spectra observed by the RGS1 
(red points) and RGS2 (blue points) with respect to Model A (RGS; see 
caption of Table[T]and Sec. l3.2K when the Gaussian emission feature cen- 
tred at 0.99 keV is removed from the model. 



tirely compatible with those found by EPIC-pn {E\ = 0.97(1) keV, 
<Ti = 0.12(2) keV, N-i = 4.7(2) x 10" 3 cm" 2 s"\ EWi = 45(1) 
eV) 

To model RGS spectra, we considered only first order spec- 
tra in 0.33-2.07 keV band, since the second order added no infor- 
mation. We modelled the continuum as in model A, with a mul- 
ticoloured disc black-body and a Comptonized component; given 
the limited bandwidth covered by the RGS, we fixed the value of 
the asymptotic power law index of the Comptonized component to 
the value indicated by the EPIC-pn analysis (see leftmost column 
of Table [T). We modelled interstellar absorption with TbNew, and 
let a relative normalisation constant between the two RGS spec- 
tra free to vary. As residuals appeared around the O I K edge at 
(0.54 keV), and close to the Ol ls-2p transition (0.5275 keV), 
we let the abundance of O as a free parameter, obtaining a best- 
fit value of Ao — 2.11 ± 0.15 with respect to a H column den- 
sity of TVh = 1.0(5) x 10 21 cm 2 . We found an excess around 
~ 1 keV in RGS data as well (see Fig. [4] where the residuals ob- 
tained removing the lines from the best-fitting models, are shown); 
by modelling it with a Gaussian feature, parameters entirely com- 
patible with those obtained modelling the EPIC-pn data with model 
A, were found. Including the Gaussian feature the chi-squared of 
the fit decreased by A% 2 = —181.2. To check the significance 
of this excess, we applied the same method of posterior predictive 
p-values described in the previous section. This yielded a proba- 
bility of p ^ 10" 3 of the observed improvement being due to 
chance. The best fitting values of the parameters are listed in the 
column labelled as Model A (RGS) of Table Q] A possible corre- 
lation with the composition of the medium yielding interstellar ab- 
sorption could be excluded, as even varying the abundances of the 
atomic species with edges falling in the considered energy band, a 
chi-squared similar to that obtained with the Gaussian feature could 
not be obtained; at the same time, contribution from absorption of 
O VIII, could be excluded as the fit with an edge fixed at 0.87 keV 
did not determine a fit of a comparable quality. 

3.3 Simultaneous RXTE - XMM-Newton spectrum 

Analysis of the observation performed by RXTE simultaneously to 
XMM-Newton allows to constrain on a broader energy band the X- 
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Figure 5. 0.33-50 keV unfolded spectrum of HETE J1900.1-2455 obtained 
fitting simultaneously spectra from RGS1/2 (0.33-2.07 keV, blue and green 
points, respectively), EPIC-pn (0.6-11 keV, red), PCA (3-30 keV, cyan) 
and HEXTE (20-50 keV, magenta) and considering Model B (see Table[TJ, 
whose component are also over-plotted (top panel). Residuals with respect 
to Model B (bottom panel). 



ray spec trum of HETE J19 00. 1-2455. By using the same bands de- 
fined bv lWatts et al. ||200i) we obtained values of the soft and hard 
colour of 1.1 and 1.0, r espectively. Compa ring these values with 
those plotted in Fig. 4 of IWatts et al] (2009), we conclude that dur- 
ing this observation HETE J1900. 1-2455 was in the softer end of 
the hard island state. 

We fitted the spectra observed by PCA (3-30 keV) and 
HEXTE (20-50 keV) on-board RXTE simultaneously to those ob- 
tained by RGS1/2 (0.33-2.07 keV) and by EPIC-pn (0.6-11 keV), 
fixing the normalisation constant of the EPIC-pn spectrum to one, 
and letting a relative normalisation constant free to vary for the 
other spectra. We obtained a reduced chi-squared of 1 .24 over 4458 
degrees of freedom by using Model B (see Table [T}, and letting 
the abundance of O in the interstellar medium as a variable pa- 
rameter. Unlike the case described in Section |3~T| we significantly 
detected a diskbb component even with a reflection component 
included in the model, as its addition decreased the chi-squared by 
A X 2 = -432 for two d.o.f. more. The Gaussian-shaped feature at 
~ 1 keV was described by parameters compatible with those found 
when modelling both EPIC-pn and RGS data with model A. The 
detection of a disc reflection component was significant also in the 
broad-band spectrum, since removing it from the model yielded 
a chi-squared increase of Ax 2 = 314 for five d.o.f. more. The 
parameters of disc reflection are compatible with those indicated 
by the analysis of EPIC-pn data alone, indicating an inclination of 
i — (30.6^2 e)° an( ^ an mner disk radius of R ln — 25^^ R 9 . 
The best-fitting parameters are shown in the rightmost columns Ta- 
ble!]] the unfolded spectrum is plotted in the upper panel of Fig.[5] 
and residuals with respect to the best-fit model shown in the lower 
panel. 

To check if a slightly mis-calibrated response of the EPIC- 
pn at low energies could affect model parameters, and in partic- 
ular those describing the reflection component, we also performed 
spectral fitting after removing EPIC-pn data below 2 keV. The best- 
fitting model had a reduced chi-squared of 1.22 (over 4427 d.o.f.) 
and no significant difference with respect to the parameters listed 
in the rightmost column of Table [T] was found. It is worth noting 
that in this case, the excess at w 1 keV was satisfactorily mod- 
elled by a reflection model with a Fe abundance of 2.0^0 20 times 



the Solar value, without the need of an additional Gaussian feature. 
At the same time, given the slight mismatch between the PCA and 
the EPIC-pn residuals (see bottom panel of Fig. |5), we removed 
the PCA data below 8 keV, but still no significant variations were 
found in the model parameters. 



4 TEMPORAL ANALYSIS 

In order to search the EPIC pn time series for a periodicity at 
the 377.3 Hz spin frequency of HETE J1900. 1-2455, we first re- 
ported the photon arrival times to the Solar system barycentre us- 
ing the_source coordinates of the optical counterpart estimated by 
lFoxl j2005). Sinc e the pulse amplitude of AMSPs usuall y decreases 
below 2 keV iPatruno et all 120091 ; IPapitto et all l201fj) . compati- 
ble with the increasing contribution of the accretion disc to the 
overall flux, we restricted our analysis to the 2-11 keV energy 
band. We corrected time delays induced by the orbital motion of 
the source in the binary system using ephemeris determined by 
IPatrunol ( l2012t) . who derived a timing solution valid over a time 
interval of 2.6 yr (MID 53539-54499) based on RXTEfPCA data. 
To ascertain whether the uncertainties on this set of parameters are 
small enough not to affect coherence of the signa l over the EPIC 
pn exp osure, we considered the relations given by C aliandro et al.l 
020120 . who estimated the reduction of the power obtained with a 
Fourier transform of a coherent signal, e 2 , as a function of the dif- 
ference between actual orbital parameters and the estimates used to 
correct photon arrival times. We evaluated the relations they de- 
rived (see their Eq. 52, 64, 68, 72) for the orbital and spin pa- 
rameters of HETE 11900.1-2455 and a signal power reduction of 
e 2 = 0.8; to produce a similar power decrease, the estimates of 
semi-major axis of the NS orbit, eccentricity, epoch of zero mean 
anomaly, and orbital period used in arrival time corrections should 
differ by Sx = 1.7 x 10" 3 lt-s, 8e = 0.04, ST* = 30 s, 
SPoib = 2.7 s, from the actual values, respectively. All but ST* 
are more than t wo orders of mag nitude larger than u ncertainties 
quoted by both lPatrunol d2012l) and lKaaret et al.l J2006l) . indicating 
how available estimates were accurate enough not to induce loss of 
signal power. The uncertainty on the epoch of zero mean anomaly 
was evaluated by propagating the uncertainty on the estima te eval- 
uated at the reference epoch of the timing solution given bv lPatrunol 
d2012l) (t = MJD 53538.76, Patruno, priv. comm.), through the 
At = 6.25 yr elap sed until the XMM-Newton observation (see, e.g. 
IPapitto et aUl2005h 
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where P or b = 4995.2630(5) s is the orbital period, o T * = 0.8 
s and op oTb = 5 x 10 -4 s are the uncertainties on the epoch of 
zero mean anomaly and orbital period, respectively, and P or h < 
1.2 x 10~ 10 is the upper li mit on the orbit al period derivative at 95 
per cent confidence level jPatrunoll2012l) . As long as such loose 
upper limit on P or t is considered, or*(tEPN) — 15 ST*, and 
nx* — 90 corrections with values of T* differing by ST* , should 
have been performed to cover all the possible value of T*, ranging 
fromT*(t E pN) -3<7t* (*epn) to T*(t E pN) + 3(tt* (*bpn). How- 
ever, t he orbital period derivative measured from SAX J 1808 .4— 
3658 dHartman et all 120081; 151 Salvo et al.l l2008t lHartman et all 
120091 ; iBurderi et alj|2009t IPatruno et alj|2012l) . a source with or- 
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bital parameters similar to those of HETE J1900. 1-2455, is ~ 
25 times smaller than the upper limit given by Patmno (2012) 
for HETE J1900. 1-2455. The value measured for SAX J1808.4- 
3658 is already so large that largely non conserv ative scenarios 
of m a ss transfer should b e invoked to explain it d5i Salvo et al.l 
2008; Burde ri et al. 1 12009b . Consider an orbital period derivative 
equal to the value measured for SAX J1808.4-3658, we obtained 
0"t*(£epn) = 23.9 s; only n' Tt — 6 <7t* (tepN)/5T* ~ 5 trials 
were then needed to cover possible values of T"*(iEPN), when a 
plausible value of the orbital period derivative of HETE J 1900.1- 
2455 was considered. 

To determine over which frequency range a meaningful search 
for a signal at the spin period of HETE J1900. 1-2455 should have 
been performed, we evaluated the spin frequency variation Av 
driven by the accretion between 2005 and 2012. Assuming the 
base level spin frequency derivative estimated by Patruno (2012), 
ii — 4.2(1) x 10~ 14 Hz s _1 , compatible with the torque ex- 
pected to be imparted by mass accretion at a rate of few x 10~ 10 
Mq yr _1 onto a NS with magnetic field ~ 10 8 G, we obtained 
Av = 8.3(2) x 10 -6 Hz. Since the frequency Fourier resolution 
of a power spectrum produced over the 60 ks EPIC pn time series is 
Svf = 1/Tobs = 1.67 x 10 Hz, one frequency bin contains both 
the value of the spin frequency in 2005 (377.29617188 Hz, Pa- 
truno, priv. comm.) and the one expected in 2012 (377.2961802(2) 
Hz). We did not detect any significant signal in power spectra of 
the time series corrected with the ri T , — 5 plausible values of 
T*, inspecting a single frequency bin centred at 377.2961735 Hz, 
of width equal to 8uf- The maximum Leahy power we obtained 
is Pmax = 4.7, lower than the 3-a confidence level threshold 
dvan der Klislll989l ; IVaughan et alJI 19941) 



P 3CT (n tot = 5) = 21n(n tot ) + 2 In = 15.0. 



(2) 



This threshold is obtained by evaluating the power level that has 
a probability lower than p = 2.7 x 10~ 3 of being exceeded by 
noise powers distributed as a % 2 with two degrees of freedom, tak- 
ing into account the five trials performed (n' T * — 5, n v — 1). We 
took into account the effect of noise-signal interactions due to the 
fact that the total power in a frequency bin is the square of the sum 
of noise and signal Fo urier amplitude, each of which is a complex 
number I Grothl 19751) . To this aim we used the routine given by 



IVaughan et alJI 19941 ; the maximum observed power translates into 
an upper limit on the signal power of Ptjl = 23.2, at 3<j confidence 
level. Considering that the total number of photons in the time se- 
ries is JV 7 = 3.96 x 10 6 , and that the averag e power reduction due 
to frequency and T* binning is = 0.773 dVaughan et al.|[l994l) 
and £t* — 1 — e 2 / 2 = 0.9, respectively , the 3-a upper limit on the 
pulse amplitude is (Vaughan et al. 199J) 
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Here, sinc(x) = sin(a;)/a;, and the Nyquist frequency is itf y 

' 1 n bin t res ) = 2117 Hz, as the time resolution of the series is 
t r es — 2.952 x 10~ 5 s and we set ribm = 16. 

Even by extending the parameter space over which we 
searched for a signal, inspecting all the values of T* compatible 
with the loose available upper limit on the orbital period deriva- 
tive (n' T , = 90, n v = 1, P 3 a{n t ot = 90) = 20.8), and the 
two nearest Fourier independent frequencies (n' T , — 90, n„ = 3, 
P3cr{ntot = 270) = 23.0), no detection w as achieved (P m ax 
1 1 . f ) ) As some of the detections claimed by iPatrunol (12012) were 
achieved on short time intervals, we performed a pulsation search 
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Figure 6. Power spectrum obtained by averaging 1940 intervals, each ~31 
s long, and re-binning the resulting spectrum with a factor 1 .03. A white 
noise level of 1.992(2) was subtracted. The solid line is the best-fit model 
composed of four Lorentzian components, which are also shown as dotted 
lines (top panel). Residuals in units of a with respect to the best-fit model 
(bottom panel). 



over n s = 120 data segments as short as 500 s (n' T * = 5, n„ = 1, 
Pza(jitot = 600) = 24.6), but still found none significant at a 3a 
confidence level (Pmax = 19.7). 

To study the aperiodic properties of the source, we averaged 
1940 power density spectra produced over ~31 s long intervals, re- 
taining photons in the full 0.3-11 keV bandwidth, and re-binning 
the resulting spectrum by a factor 1.03. Four Lorentzians compo- 
nents were used to model the power spectrum obtained after the 
subtraction of a white noise level of 1.992(2) Hz" 1 (see Fig.©: 
two flat-top noise components (i.e. centred at zero frequency), with 
width Wi = 0.80(1) and W 2 = 47(7) Hz, and two narrower 
Lorentzians centred at v 3 = 1.98(2) and va = 2.6(2) Hz, with 
width W 3 = 0.4(1) and W A = 3.6(3) Hz, labelled as L u L 2 , 
Li and L4, respectively. A similar decomposition is remarkably 
simila r to that shown by at oll sources in the hard st ate dOlive et al.l 
Il998l) . and other AMSPs dvan Straaten et al]|2005t) . In particular, 
the frequency and quality factor (defined as v/W) of features we 
denoted as L3 and L4 are similar to those of the quasi periodical os- 
cillations identified bv lvan Straaten et al I d2005l) as Llf and Lh, in 
low luminosity, non-pulsing atolls, such as IE 1724-3045 and GS 
1826-24. High frequency quasi periodical oscillations were also 
searched for, by averaging 16 s long intervals, but none was found. 



5 X-RAY BURST 

A type-I X-ray burst was observed during XMM-Newton exposure 
(see inset of Fig.Q]for the 0.6-1 1 keV light curve), starting on MJD 
55823.69344. During first ~ 10 s since the burst onset, EPIC-pn 
count-rate saturated at a level of ~ 1900 s _1 due to telemetry lim- 
itation, and the light curve cannot be considered representative of 
the true burst light curve shape (see top panel of Fig. [7). Presum- 
ably for the same reason, data were not acquired during intervals 
14 - 24 and 34 - 44 s, since the burst start. At later times since 
the burst onset, the profile could be modelled with an exponential 
decay with a time constant of r = 46(1) s. This is significantly 
longer than the length of the majority of bursts observed by RXTE 



and H ETE dSuzuki et al]|2007l ; ballowav et alj|2008l : 
2009), except for the burst labelled as B 1 in Fig. 2 of 



Watts et al 



Watts et al 
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Figure 7. 0.6-11 keV light curve of the burst observed by EPIC-pn on- 
board XMM-Newton. Time is shown in seconds since the burst onset. 
The horizontal bar indicates the time interval during which a signal is 
barely detected at a Leahy power of 21.6. The inset shows the profile ob- 
tained by folding data obtained during such interval, around the frequency, 
377.04(25) Hz. 



l2009t and which also took place where the source was in the hard 
state. However, there does not seem to be a clear correlation of burst 
duration with the position of the source in a colour-colour diagram; 
the position of the source during the observation considered here 
is in fact much closer to that sh own by the sour ce before shorter 
bursts (e.g. burst B4 in Fig. 4 of lWatts et al.l20 09). than at the onset 
of the other long burst detected. 

We searched for burst oscillations over 4 s intervals of EPIC- 
pn data; we inspected frequencies ranging from 374 to 378 Hz, for a 
total of 16 trials in every interval, in order to cover the i nterval over 
which the frequency of burst oscillations discovered bv lWatts et al.l 
( 2009) was observed to vary. We found a barely significant signal 
at a frequency of 377.04(25) Hz with a Leahy power of 21.6, in the 
interval starting 24 s after the burst onset. Such a power is larger 
than the 3<r threshold for a single interval (17.2), but reduces to a 
significance of 2. la if the number of intervals considered to analyse 
the burst (100) is taken into account; we cannot therefore exclude 
it represents a statistical fluctuation. The folded profile is shown in 
the inset of Fig. [7] If the oscillations were effectively present, this 
would be the first time t hey were observed in the hard state, since 
oscillations reported by IWatts et all ((2009) were found in a short 
duration burst, ignited when the source was in a soft state. We note 
that before the interval where the marginal detection was achieved, 
sensitivity to pulsed signals was severely reduced due to telemetry 
overload. 

For the same reason, time-resolved spectroscopy of the burst 
could be meaningfully performed only after ~ 25 s since the burst 
onset. No background was subtracted, and the best-fit Model A (see 
leftmost column of Table[TJ was used to model persistent emission, 
letting only the absoiption column free to vary. Spectra in the 0.6- 
1 1 keV energy band were satisfactorily fit by an absorbed black- 
body, with a temperature decreasing from 1.6(1) to 0.9(1) keV at 
the end of the burst, and an apparent radius roughly constant around 
a value of ~ 6 ds km. The value of the absorption column density 
we found, iV H = (0.99 ± 0.15) x 10 21 ctrT 2 is compatible with 
estimates indicated by analysis of persistent emission (see Table[T](. 



6 DISCUSSION AND CONCLUSION 

LMXB hosting a NS are usually classified in Z and atoll sources, 
based on the track follo wed in a colour-colour diagram and the ape- 
riodic time variability (Hasing er & van der Klislll989), with atolls 
emitting a luminosity ,$ 0.5 Lbm, and showing a harder spec- 
trum, on average. Similar to other AMSP, HETE J1900. 1-2455 be- 
haves as a typic al atoll source, sp ending most of the time in the 
hard island state (Watts et al. 2009). Also the observations consid- 
ered here caught the source in such state, while emitting a 0.33-50 
keV unabsorbed luminosity of 4.07(3) x 10 36 erg s _1 ; roughly 
90 per cent of the X-ray observed in such band is described by a 
F ~ 2 hard component, with no cut-off significantly detected in 
our data. To physically model such a component we used Comp- 
tonization of soft, kT so [ t ~ 0.45 keV, photons in a medium with 
an electron temperature fixed to 100 keV. The slope of such hard 
component and the a ssumed electron temperature indicate an op- 
tical depth, t — 1 dZdziarski et alj|l996t Evcki et all 1 19991). A 
similar Comptonized component was used bv lFalanga et al.l l 20071) 
to model the 2-300 keV spectrum observed by RXTE and INTE- 
GRAL during observations performed in 2005, and is invariabl y 
found to dominate the X-ray spectrum of AMSP (Poutanen l2006t) . 
Considering also the e nergy-dependent properties of pulse profiles, 
iGierlinski etail d2002l) interpreted this component as originated in 
the accretion columns, where mass is channelled by the magnetic 
field and heated by the interaction with radiation coming from the 
NS surface. However, the non detection of pulsations casts doubts 
on the presence of accretion columns during the observation of 
HETE J 1900. 1-2455 presented here. 

Two broad emission features centred at energies 0.98(1) and 
6.58(7) keV were detected in EPIC-pn spectrum. The energy of 
the latter is com patible with Kq trans ition of Fe XXIII-XXV (6.6- 
6.7 keV, see, e.g., Kal lman et alj2004Q; such detecti on confirms the 
results previously obtained by Cackett et al. l l2010l) with a Suzaku 
observation, and by iFalanga et al.l 1 120071) , even if at a much lower 
spectral resolution. Origin of the ~ 1 keV feature is instead less 
clear. Calibration uncertainties of the EPIC-pn spectra taken in tim- 
ing mode were often invoked to explain appearance of strong fea- 
tures at energies £ 1 keV; however, the presence of a similar line- 
shaped residual also in spectra observed by EPIC-MOS2 and RGS 
strongly suggests a physical origin. Both the features have a broad- 
ness a/E ~ 0.1, compatible with a similar mechanism determin- 
ing their shape. The easiest interpretation of such width is in terms 
of reflection of hard X-ray photons in the accretion disc, where 
the strong velocity field broadens discrete features, and relativistic 
effects distort and smear their shape, shifting them towards lower 
energies. In the case of the observations analysed here, a mildly 
ionised £ ~ 750 erg cm s _1 disc reflection component, broad- 
ened by motion in an accretion disc and representing 10 per cent 
of the overall flux, successfully models the iron line and improves 
description of the spectral continuum. However, a self-consistent 
modelling of the two discrete features, could not be obtained, as far 
as the EPIC-pn data below 2 keV are considered. Strong emission 
features at energies similar than those found here were recently re- 
ported thanks to XMM-Newton spectra of a couple of narrow-line 
Seyfert 1 galaxies, and interpreted in terms of emission lines due 
to Fe Ka and Fe La transitions takin g place as the accretion disc 
reflects hard illuminating emission (Fabian etal1l2009l; I Ponti et al.l 



l201fJ ; IZoghbi et alJl20ld;lFabian et alj|2012l) . Despite the ratio be 
tween the photon fluxes in the two features observed by EPIC-pn 
(= 6.8 ± 1.7) is less extreme than those they reported ( ^ 20), it is 
still too large to be described by a single reflection environment like 
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that we considered. A lower flux ratio is in fact expected when the 
illuminating spectrum becomes softer than T ;$ 2 at a given flux, 
and iron is less ionised dRoss & Fabiarj2005h , However, excluding 
from the fit EPIC-pn data below 2 keV and modelling low energies 
with RGS data alone, we could find a self-consistent modelling of 
the two features in the same reflecting environment, with an iron 
abundance ~ 2 times larger than the Solar value. Another, not nec- 
essarily alternative interpretation of the excess is in terms of Ka 
transition of overabundant helium-like Ne; though, it was not pos- 
sible to vary the abundance of this element in the reflection model, 
and simply adding a blurred feature at the relevant energy did not 
result in a completely satisfactory description. The failure in deriv- 
ing a simultaneous description of the feature detected by EPIC-pn 
possibly follows from superimposition of different transitions not 
completely taken into account by the reflection model used, differ- 
ent abundances of ionised elements, or because of residual calibra- 
tion uncertainties. The satisfactory modelling of the two features 
by a reflection component with overabundant iron, obtained when 
EPIC-pn data at low energies are discarded, favours the latter hy- 
pothesis. In any case, the simultaneous detection of such feature 
from different instruments makes it a very promising additional 
probe of the accretion flow to LMXB hosting NS, the nature of 
which will be firmly assessed by future observations. 

Fitting the smeared refle cted component with the rdblur 
kernel of iFabian et al.l dl989h . we estimated the inclination, the 
emissivity, the inner radius of the disc, i — (30. 61?, P — 
— 3.7l2'g> ^ in — 25l^ R g , respectively. From Doppler tomog- 
raphy of a broad Ha line assumed to originate in the accretion disc, 
lElebert et al.l d2008l) estimated the inclination of the system to be 
i 20° for a 1.4 M©, a value not far from that reported here. Our 
estimate of the inn er disc radius i s co mpatible with the value of 
14(1) R 9 found bv lCackettet"al] ( l2010L for a rather low inclina- 
tion of 4°), who modelled a Suzaku spectrum of this source taken 
when its 0.5-25 keV luminosity was ~ 5 x 10 35 erg s _1 (roughly 
a factor three times lower than during the observations presented 
here) by adding a reflection component to a T ~ 2.2 power law. At 
a distance of ~ 25 R g from the compact object, relativistic effects 
are small and the line shape we observed is rather symmetric. How- 
ever, alternative interpretations proposed to explain broad symmet- 
ric iron lines, such as Compton broadening in a hot corona, are not 
compatible with what observed in the broadband spectral decom- 
position. An optically thin, r ~ 1 corona cau ses an average photon 
energ y shift of Ae/e = (4fcT e — e)/m e c 2 JRvbicki & Lightmanl 
1979); to produce the observed line width of 0.7(1) keV, the elec- 
trons in the Comptonizing cloud should have a temperature of ~ 15 
keV and a similar component was not observed in the spectrum. 
The value estimated for the inner disc radius translates into a value 
of 52^23 km, for a 1.4 M NS. Despite the large uncertainty af- 
fecting this estimate (Ri n > 12 km at 99 per cent confidence level), 
the shape of reflected spectrum indicates that the accretion disc is 
truncated. A similar size of the inner rim of the disc is suggested 
also by the kT[ n ~ 0.2 keV thermal component detected at low 
energies; if modelled with an accretion disc model, its normali- 
sation indicates an apparent inner disc radius ~ 30 ds km (for 
i = 30°). The actual inner radius of the disc can be larger up to 
a factor of two than this estim ate, when colour corrections are con- 
sidered dMerloni et alj :2000) and torque is assumed to be retained 
through the tr ansition region from th e optically thick disc to the 
inner hot flow ( lGierlinskietailll999h . At the same time, consider- 
ing the luminosity of the photo-ionising Comptonized component 
(Lirr ~ 3.5 x 10 36 d 2 erg s _1 ), the ionisation parameter of the 
observed reflected component (£ = Li rr /nnp 2 — 750 erg cm s _1 ) 



is expected to be produced at a distance p ~ 30 km from the source 
of illuminating photons, for a hydrogen density of 5 x 10 20 cm -3 , 
a value appropri ate for a Shakura S unyaev disc around an accreting 
NS (see, e.g., iRoss & Fabianll2007l) . Further, for a 10 km NS and a 
system inclination of 30°, the observed ratio between the reflected 
and illuminating component (~ 0.1, corresponding to a reflection 
amplitude Q./2tt ~ 0.13) ind icates an inner disc radius £ 50 km 
( Gierli hski & PoutanerJl2005l who calculated the reflection ampli- 
tude as a function of the inner disc radius in the case of a small 
illuminating spot at the rotational pole of the NS; a more extended 
source of illuminating photons produces a larger reflection ampli- 
tude at a fixed disc radius). 

The main peculiarity of HETE 11900.1-2455 is coherent pul- 
sations intermitt ence, whose causes are not yet fully understood. 
ICummingl |2008) proposed that intermittence is caused by a screen- 
ing of the magnetosphere and burial of the magnetic field under 
the NS surface. O n the other hand, iRomanova et al] d2008l) and 
lLamb et alj d2009l) . while involving the presence of a magneto- 
sphere, explained intermittence in terms of penetration of in-falling 
matter through the magnetic field lines, and of movements of a hot 
spot located close to the spin axis of the NS, respectively, depend- 
ing on variations of the mass accretion rate. Alternative scenarios 
proposed to explain the paucity of pulsars among LMXB, and pos- 
sibly related to pulse intermittence, include scatteri ng in a hot, rel- 
ativel y optically thick cloud surrounding the NS dTitarchuk et al] 
2002), a nd red uction of pulse amplitude due to gravitational light 
bending dOzell 12009). However, considering how a Comptonized 
component produced in an optically thin, r ~ 1, corona was ob- 
served regardless of pulse detection, and how the compactness of 
the NS hardly varies on timescales of few years, these two lat- 
ter interpretations do not reproduce easily the behaviour observed 
from HETE 11900.1-2455. The evidences we obtained from spec- 
tral analysis of an accretion disc truncated quite far from the NS 
have to be discussed in view of the time variability shown by the 
source. No pulsations were detected during the persistent emis- 
sion, with a very low upper limit of 0.4 per cent, at 3 a confi- 
dence level, on the 2-11 keV pulsed fractional amplitude. Such 
a limit i s of th e same order of the weaker pulsations detected b y 
IPatrund d2012h . As already suggested by iGallowav et all d2008l) . 
this source could persistently show pulsations at a very low, ,$,0.1 
per cent, amplitude, beyond the detection sensitivity of current in- 
struments. The inner disc radius derived by spectral modelling of 
the reflection component is compatible with models in which a 
magnetosphere is able to truncate the optically thick disc, and co- 
herent oscillations are reduced in amplitude by geometrical effects 
dLamb et al1l2009l) and/or b y wea k channelling of accreting mat- 
ter, if any I Romanova et al] [20081) . However, the tight upper limit 
on pulse amplitude and the absence of a clear relation between ap- 
pearance of pulses and variations of the mass accretion rate, do not 
necessarily favour them. On the other hand, spectral and aperiodic 
timing properties of this source are entirely indistinguishable to that 
of (assumed) non pulsing atolls in the hard state, suggesting how 
the accretion flow to HETE 11900.1-2455 is probably not chan- 
nelled by a magnetospher e. This would ke ep open the possibility 
of a buried magnetic field dCummin J2008h . In such a scenario, the 
optically thick flow in a geometrically thin accretion disc would be 
replaced by a inner optically thin hot flow, rather than by a magne- 
tospheric flow. This fits into the truncated disc scenario explaining 
hard state s shown by both accreting black holes and neutron stars 
(see, e.g. jDone et al.|[2~007l . for a review). Similar estimates of the 
inner disc radius were indeed obtained from iron line modelling of 
non-pulsing atolls in the hard state (e.g., the cases of 4U 1705-44 
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and MXB 1728-34; IeTAi et alj|20ld : lEgron et alj[201 ll l2012h . It 
remains an open question whether the accretion disc flow is trun- 
cated by a similar mechanism, connected to the overall evolution of 
the accretion disc rather than to the NS magnetic field, even when 
a magnetosphere is certainly present (i.e., when pulsations are de- 
tected). Long observations of AMSPs performed while they show 
pulsations, and aimed at a detailed modelling of broad emission 
features, will be crucial to this end. 

During the XMM-Newton observation a type-I X-ray burst was 
detected. Because of telemetry overflow, we could not determine 
the peak burst luminosity, nor if a photospheric radius expansion 
took place. We barely detected a signal at a frequency of 377.04(25) 
Hz during a 4 s interval, 24 s after the burst onset; however, the ob- 
served power is low, and the signal is significant only if the number 
of burst intervals inspected is not taken into ac count. Burs t oscil- 
lations from this source were already found by IWatts et all d2009t) 
during a short burst (r = 7-8 s). The burst during which they 
detected pulsations ignited when the source was in the soft (ba- 
nana) state, with a a 2-16 keV flux of « 63 mCrab as observed 
by RXTEfPCA, which is to date the highest persistent flux at which 
bursts have been detected in this source. Oscillations were detected 
during the initial stages of the burst decay, at a frequency ~ 1 Hz 
below the spin frequency of the source, drifting towards this value 
as the burst evolved. On the other hand, the 46 s time decay of the 
burst reported here, the persistent source flux (24 mCrab in the 2- 
16 keV range) and the hard spectral state are not easily reconciled 
with the general picture of the non-pulsing LMXBs and intermit- 
tent pulsars, which preferably show oscillations during short bursts 
ignited when the source is i n the soft state, pres umably accreting 
at a rate larger than average (IGallowav et al . 2008). Moreover, the 
initial burst oscillation frequency would be within few tenths of 
Hz from the spin frequency of the source, 377.296 Hz, as it is 
commonly observed for oscillations occurrin g during the decay of 
bursts shown by persistently pulsing AMSP (Watts 2012); though, 
as the oscillations barely detected here were observed at later times 
since the burst onset, such frequenc y could be compati ble with the 
positive frequency drift observed bv lWatts et al.l d2009l) . Assuming 
burst oscillations marginally detected during this observation to be 
real, they would be much more similar to those of persistently puls- 
ing AMSP, in terms of spectral stat e at the burst ignit ion and fre- 
quency drift, than those detected bv lWatts etaUT 2009); this could 
indicate a role played by a residual magnetic field during the burst 
reported here, under the hypothesis that the still poorly understood 
properties of burst oscillations are related to the channelling of ac- 
creted material by the NS magnetic field (see, e.g., the review by 
IWattsll2012h . 
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